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Abstract 

We generalize the computation of supersymmetric contributions to the muon anoma- 
lous magnetic moment (g — 2)^ to the NMSSM. In the presence of a light CP-odd Higgs 
scalar, these can differ considerably from the MSSM. We discuss the amount of these 
contributions in regions of the parameter space of the general NMSSM compatible with 
constraints from B physics. In the mSUGRA-like cNMSSM, constraints from [g — 2) M 
prefer regions in parameter space corresponding to a low SUSY breaking scale. 
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1 Introduction 



The result of the measurement of the anomalous magnetic moment of the muon = (g^ — 2)/2 
by the E821 experiment at the Brookhaven National Laboratory (BNL) [1] can be considered 
as a possible hint for physics beyond the Standard Model (SM). The determination of the 
SM contributions to requires - amongst others - to compute the leading order hadronic 
contribution a^ LO , which is the main source of the present SM uncertainty. Recent calculations 
of a^ LO via hadronic e + e~ data [2-6] are in good agreement, and lead to ~ 3 standard deviations 
between the experimental value for a M and the SM prediction. 

Alternatively, a^ LO can be determined from hadronic r-decays [7,8], leading to an agreement 
within one standard deviation between the experiment and the SM prediction for a^. However, 
a comparison of the r branching fractions into pions with the corresponding e + e~ spectral 
functions reveals a discrepancy of ~ 4.5 standard deviations [6] and requires assumptions on 
the pion form factor, isospin violating effects and vector meson mixings. Since e + e~ data is more 
directly related to a^ LO , it is advocated to use the hadronic e + e~ data only for a computation 
of al LO [5]. 

Additional contributions to a M appear within supersymmetric extensions of the SM (see [9] 
for a recent review). One-loop [10-13] and two-loop [14-22] contributions have been evaluated 
within the Minimal Supersymmetric Extension of the Standard Model (MSSM). The conclusion 
is that the MSSM is able to explain the 3 standard deviations between the experiment and the 
SM, provided that the masses of the electroweakly interacting supersymmetric particles are 
not far above the electroweak scale [9, 12, 13, 19, 20]. On the other hand, the measured value 
for a M provides constraints on the parameter space of the MSSM such as the positivity of the 
supersymmetric Higgs mass parameter \i. 

The purpose of the present paper is the study of constraints from the measured value for 
a M on the parameter space of the Next-to-Minimal Supersymmetric Extension of the Standard 
Model (NMSSM). First estimates of the contributions to a M in the NMSSM have been performed 
in [23] (see also [24] for a M in a similar U(l)' model), but subsequently we aim at an accuracy 
comparable to the one in the MSSM (limiting ourselves, however, to dominant logarithms 
at two-loop order). To this end, most of the corresponding formulas can be translated in a 
straightforward way from the MSSM to the NMSSM. 

Nevertheless, the numerical results in the NMSSM can differ considerably from the ones 
in the MSSM: the NMSSM contains an additional gauge singlet superfield S, the vacuum 
expectation value of which generates the Higgs mass parameter fi of the MSSM [25]. The CP- 
even and CP-odd components of S will generally mix with the neutral components of the two 
(MSSM) Higgs doublets H u and Hd- Depending on the parameters of the NMSSM, the lightest 
neutral CP-odd Higgs scalar can be quite light [26] and lead to numerous new phenomena [27]. 
Such a light neutral Higgs scalar can also have an important impact on [28]. In section 3 
we will study possible effects of a light neutral CP-odd Higgs scalar on in regions of the 
parameter space of the NMSSM, which satisfy present bounds from LEP [29] and 5-physics [30]. 

In the remaining part of the introduction we briefly review the various SM contributions to 
a M in order to clarify, which additional contribution would be desirable. The SM contributions 
to a At can be split into pure QED contributions a® ED (known to four- loop order), leading order 
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hadronic contributions a BLO , next-to- leading order hadronic contributions a^ NLO (vacuum 
polarisation diagrams only) and a LBL (light-by-light contributions only), and electroweak effects 
a EW _ rj^g gum j g ^ Q De com pared with the experimental result [1] 



a EXP = 11 659 208.0(5.4)(3.3) x 10" 10 . (1.1) 
The QED contribution is [31] 

a QED = ii 6 5 8 47i. 8 n3(i 6 2) x 10~ 10 (1.2) 

(taking into account estimated five-loop contributions [32]) with an error well below the exper- 
imental one. The remaining difference is 

4 XP ~ 4 ED = (736.2 ± 6.3) x 1(T 10 . (1.3) 

Among the recent evaluations of the leading order hadronic contributions a BLO [2-6] (including 
the most recent data from SND, CMD-2 and BaBar) we use - in order to remain conservative 
- the largest estimate from [2,5], leading to the smallest deviation from the SM: 

a BLO = (692.1 ±5.6) x 1(T 10 (1.4) 

For a BNLO we use the same reference [2,5]: 

a HNLO = (_iQ.03 ± 0.22) x 10" 10 (1.5) 

The most recent determination of the hadronic light-by-light contribution [33] gives 

a LBL = (11.0 ±4.0) x 1(T 10 . (1.6) 

Adding all errors quadratically, one obtains 

< XP - a® ED+HAD = (43.1 ± 9.3) x lO" 10 ; (1.7) 

the remaining discrepancy should be explained by electroweak effects and/or contributions 
beyond the SM as supersymmetry. Within the SM, the electroweak contributions are to two- 
loop order [34] 

a EW = (15.4 ±0.2) x 10~ 10 , (1.8) 

which leads to the present discrepancy of about three standard deviations between the experi- 
ment and the SM: 

4 XP - a T = ( 27 - 7 ± 9 - 3 ) x 1(rl ° (1.9) 
In the next chapter 2 we will review the supersymmetric contributions to a M , and specify 
their dependency on the parameters of the NMSSM. Contributions to a M which are very similar 
in the NMSSM and the MSSM (without possibly light CP-odd Higgs scalars), and contributions 
which involve the possibly light pseudoscalar of the NMSSM, are treated in separate subsections. 
The formulas of chapter 2 will be made public in the form of a Fortran code on the NMSSMTools 
web page [35]. In chapter 3 we present numerical results for a M in various regions of the 
parameter space of the NMSSM, amongst others regions with a light CP-odd neutral Higgs 
scalar in the spectrum, and in the cNMSSM (with universal soft terms at the GUT scale). At 
the end of chapter 3 we conclude with a short summary. 
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2 Supersymmetric Contributions to in the NMSSM 



2.1 Contributions without pseudoscalars 

The one-loop contributions to a M in supersymmetric models are known to consist in chargino/ 
sneutrino or neutralino/smuon loops [10-13]. The corresponding expressions are the same in 
the MSSM and in the NMSSM, provided that the additional singlino state in the neutralino 
sector (which includes now five states) is taken into account. Using the formulas in [11], the 
chargino/sneutrino and neutralino/smuon contributions to a M can be written as follows: 
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[2.2) 



Here , c£, nf m and n\ m denote the couplings of the mass eigenstates of the charginos 
or neutralinos to the sneutrino or smuons and the muon, which depend on the corresponding 
mixing matrices described below: 
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(2.3) 
(2.4) 
(2.5) 

(2.6) 



The conventions for the mass matrices and the resulting mixing matrices are as in [11,35-37]: 

• The neutralino mass eigenstates x^i i = 1, . . . , 5 in the NMSSM (ordered in increasing ab- 
solute mass) are given in terms of the interaction eigenstates ipj = (— z'Ai, — iX 2 , ^>°, if) 8 ) 
by \ { ! = Nirfl 

• The chargino mass eigenstates x k i ^ = are related to the charged gaugino and 
higgsino interaction eigenstates ip + = (— zA + ,?/>+) and ip~ = through the 
rotation matrices U, V: xt = Vki^t ', Xk = Uki^pf ■ 



The smuon mass eigenstates jl m , m = 1, 2, with masses m^ m , result from the interaction 
eigenstates \x l = (Jxl, Jxr) and the rotation jl m = X^jl^. In the following, this definition 
of the matrix X% m will be extended to all sfermions /, in which case it will be denoted 
as XL. 



The /z-sneutrino mass is written as mj 



x ; 



and Xk denote the following mass ratios: 



m 2 n/m 
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Am ' 



m 



jm\ . The functions F 1 ! l N are given in the appendix, is the 
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muon Yukawa coupling: — ^p, where is the vacuum expectation value of the Higgs 
doublet which couples to down quarks and leptons. 

Since contributions to a M require a chirality flip, the formulae (12.11 12. 2p involve terms - 
apart from the prefactor m M - which are proportional either to the muon mass (if the chirality 
flip occurs in the external legs) or to a chargino/neutralino mass (when it occurs in internal 
lines). In practice, the numerically dominant contributions usually originate from a chirality 
flip in the internal lines of the chargino/sneutrino contribution, i.e. the second term in (12.11) . 
It is proportional to , and hence to the Yukawa coupling ~ m M / cos (3 ~ m M tan (3 (for 
tan/3 > 1). 

An analysis of the chargino/sneutrino diagram reveals (see, e.g., [9]) that this contribution 
carries the sign of the Higgs mass parameter [i (= fi e ff in the NMSSM), hence a positive value 
for n is phenomenologically favoured. The contribution decreases with increasing sneutrino 
mass. For positive /i, the chargino/sneutrino contribution can well resolve the discrepancy 
(11.91) . provided the muon sneutrino is relatively light or tan/3 is large; see, e.g., Figs. 1-3 in 
section 3. (For a light sneutrino and large tan/3, this contribution can even be larger than 
desired.) 

It has been pointed out in [11,12] that the neutralino/smuon contribution (12.21) could also 
explain the 3a discrepancy (II. 9p . if the bino is quite light and the smuon mass eigenstates are 
not too heavy. This bino/smuon contribution has the interesting property to be quite insensitive 
to tan (3. 

In the NMSSM, a light neutralino could also be dominantly singlino-like. However, in this 
case its contribution to a M is suppressed because of the weak couplings of the singlino to the 
MSSM sector. 

In addition to these one-loop diagrams, two-loop contributions to in the MSSM have 
been studied in [14-22]. First, we include large logarithms arising from QED corrections to 
one-loop diagrams as computed in [14]: 



This leads to a reduction by a few percents of the LO contributions. 

Additional bosonic electroweak two- loop diagrams were computed in [20], which can be 
written as 



(Up to the more complicated Higgs sector of SUSY models, (12.81) contains two-loop electroweak 
SM contributions included in (11.81) . We took care not to count the SM contribution twice.) 
Subsequently we confine ourselves to leading logarithmic contribution ~ c 2 I LBos , which reads [20] 




(2.7) 




(2.8) 
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In the MSSM, the Higgs contribution is of the form (see eq. (27) in [20]) 
h cos 2/3 Mf cos a cos(a + (3) sin a sin (a + (3) 



(2.10) 
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where and a are the CP-even Higgs masses and the mixing angle in the MSSM. In order 
to generalize the Higgs contribution to the NMSSM, it is convenient to re-write eq. (12.101) in 
terms of elements of the (inverse) CP-even Higgs mass matrix M<7 2 , here in the basis (H u , H^): 



c\ = cos2/3Mi [{Mf) 22 - tan/? (M^J 



(2.11) 



(Now it is straightforward to verify the sum rule d l L = \ [20] with the help of the tree level mass 
matrix M|.) 

Eq. (12. lip can be interpreted as the result of the evaluation of the Higgs-dependent two- 
loop diagrams, where Feynman rules in the interaction basis (H u , Hd) have been used and 
mass insertions were treated perturbatively - this procedure reproduces the leading logarithms. 
Furthermore, in the interaction basis it is easy to identify the additional contributions from the 
singlet scalar of the NMSSM, which decouples from the muon as well as from gauge bosons. 
The only possible additional contributions arise from the coupling of the singlet to charged 
Higgs bosons (involving always at least one power of the Higgs singlet-doublet coupling A), and 
the simultaneous presence of a mass insertion transforming a doublet into a singlet (also of 
0{\)). Since the coefficient of c£ in ( 12.81) is only ~ — 2 x 10~ n , we neglect subsequently effects 
of 0(X 2 /g 2 ) in c h L . Then eq. ( 12. lip is equally valid in the NMSSM; all that remains to be done 
is to replace M| by the 3x3 mass matrix of the CP-even sector of the NMSSM. Returning, 
for convenience, to mass eigenstates, c h L is then of the form 



Here, the conventions for the Higgs mixing matrices are as follows [35] (for convenience we 
discuss the complete Higgs sector here, which is useful for what follows below; note that H u 
and Hd are exchanged w.r.t. the conventions in [36,37]): 

• The CP even Higgs eigenstate hi, i = 1,2,3 (ordered in increasing mass) in the NMSSM 
are a mixture of the real parts of the neutral Higgs components S 1 = (H u r, HdR, Sr): 



• The CP odd Higgs states a^, % — 1, 2 in the NMSSM originate from the imaginary parts 
of the neutral Higgs components P 1 = (H uI , H d i, Si) (after omission of the Goldstone 
boson G° = — sin f3H uI + cos (3H d i): en = P' iX (cos [3H uI + sin (3H dl ) + P^Si- 

• The charged Higgs boson H ± is obtained from the charged Higgs components H^, 
as E ± = cospH± + sinpH^. 

Two-loop diagrams involving closed sfermion and chargino loops were studied in [18, 19] 
and [20], respectively. The leading contributions (photonic Barr-Zee diagrams) can be found in 
[9,18] in the context of the MSSM and can be generalized to the NMSSM through a replacement 
of the corresponding couplings. 




(2.12) 
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For the sfermionic diagrams we use the formula (5) of [18]: 

3 

1^ N cQ) ~2 J f 



f «=1 / 



m ■ 



m 



(2.13) 



hi 



where N[ corresponds to the number of colors (3 for squarks, 1 for sleptons), Qj is the electric 
charge, and the Higgs/sfermion couplings in the NMSSM are given by: 
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(2.16) 



while the muon/Higgs couplings are simply A^ = v u , v d and s denote the vacuum expec- 
tation values of H u , H d and S, respectively, in the normalisation where v\ + v\ = 1/(2\^2Gf)- 
A is the Higgs singlet-doublet Yukawa coupling, and A t ^ T are the soft supersymmetry breaking 
Higgs-sfermion trilinear couplings for the third generation. 

For the closed chargino loop, we employ the formula (63) of [9]: 



Gpfn 2 a 



k=l 



Re[X h u i X h 




Re[A^A ai ± ] fm 2 ± 
— J PS . 

mi 



(2.17) 
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with the chargino/Higgs couplings of the NMSSM: 

A\ = ^^[\U k2 V k2 S i3 + g 2 (U kl V k2 S a + U k2 V kl S i2 )} (2.18) 

X fc g 2 

= — — — — ^- [\U k2 V k2 P^ 2 — g 2 (U k iV k2 cos/3 + U k2 V k i sin/3) P^] (2.19) 



x 



.92 



The functions fj, f$ and /p^ can be found in the appendix. 

Among the missing contributions are one-loop diagrams with the exchange of a Higgs boson 
between the muons, and two-loop diagrams with a closed SM fermion loop involving Higgs 
bosons (and a photon). Since these can play a particular role in the NMSSM, they will be 
treated separately in the next subsection. 



2.2 Contributions including pseudoscalars 

Higgs effects are usually negligibly small in the SM or the MSSM because of the existing lower 
bounds on the Higgs masses. The SM one- loop Higgs/muon diagram is, indeed, about four 
orders of magnitude below the sensitivity of the BNL experiment for a SM Higgs mass above 
114 GeV. However, this is not necessarily the case in the NMSSM: while the lightest CP-even 
Higgs boson mass cannot be far below the LEP bound unless it decouples from the SM sector, 
the lightest CP-odd boson a± can be as light as a few GeV [26,27]. Bounds from B-physics [30], 
especially from BR(B S — > fi + fi~), can still be satisfied for low values of tan/5 or when the 
loop-induced b — s — a\ coupling is suppressed. 

In the context of two-Higgs-doublet-models, the impact of light (pseudo-) scalars on a M was 
already pointed out in [15,16,28]. In the NMSSM, a short analysis of a u has been performed 
in [23]. 

In the following we include contributions from the Higgs sector to a u up to the two-loop 
level. One-loop scalar /fermion diagrams have been known for quite a while; see, e.g., [38]. For 
completeness we detail all the SUSY Higgs contributions (CP-odd, even and charged): 

Rn \L CP even _ ^ ^ m fi \^ Sj2 C X 2 (2~x)d,X 



\L CP odd 



4^** cos' /U (1-x) 
Gujnl .„ „ „ f 1 x 3 dx 



f^E^fW/ 3 / (2.21) 

4V27r 2 tr Jo x 2 +(^) (1-x) 



G V m l 2 a f 1 X ( X - !) dx 



Sa lL charged = JLl^^ ^ ' — _ (2.22) 



For the two-loop Higgs diagrams involving a closed SM fermion loop (in contrast to closed 
sfermion/chargino loops considered in the previous subsection), we follow the analysis of [16] 
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and generalize it to the NMSSM: 



8a 



2L CP even 



5a 



2L CPc 



4V27T 3 ^ 



4 SnSi2 



3 sin /3 cos /3 



rn 



h, 



+ 



1 g| 

3 cos 2 (3 



fs 



in j 



C2 
D i2 

cos 2 /3 



fs 



3 fpS 



m. 



mi 



f tan^ If p S 



mi 



mi 



+ f 



PS 



mi 



m 



mi 



mi 



(2.23) 



(2.24) 



where the functions fs, fps are defined in the appendix. 

As noticed in [28] , one- loop and two- loop light Higgs contributions are of opposite signs and 
interfere, therefore, destructively. In the case of a CP-odd scalar, the one-loop contribution 
is negative and worsens the discrepancy (11.91) correspondingly. However, for a light CP-odd 
Higgs heavier than ~ 3 GeV, the positive two- loop contribution is numerically more important. 
The sum of both contributions is maximal around m ai ~ 6 GeV, though fairly constant in the 
range 4 — 10 GeV (see Fig. 5 below). Both one- and two-loop contributions are proportional 
to the product of two muon Yukawa couplings, which leads to an enhancement quadratic in 
tan/3. They are also proportional to the square of the mixing of the light pseudoscalar to 
the doublet sector (~ (P^ sin/3) 2 , if we consider the dominant component only). In spite 
of the appearance of P[ x , this coupling can be large enough to allow the light pseudoscalar 
contribution to reach the experimental la range of a M by itself, provided that tan f3 > 30. 
Hence, this contribution can alleviate the upper bound on slepton masses, which can be derived 
under the assumption that the chargino contribution (12. II) explains the deviation (11. 9p . 

Finally, we estimate the theoretical uncertainty following the analysis of [9], allowing for 
a 2% error on the one-loop contributions and a relative error of 30% for the two-loop results. 
Except for the light pseudoscalar contributions, which we include in the error computation 
described above, we do not expect large sources of uncertainties different from the MSSM. 
Therefore, we use the same additional constant terms as [9] for the evaluation of the error. 



3 Results 

In this section we discuss a few phenomenological aspects of the supersymmetric contributions 
to a M discussed above. First we consider the tan (3 and slepton mass dependences in NMSSM 
scenarios without a light CP-odd scalar; as expected, these results are essentially the same as 
in the MSSM. Then we examine the case of a light pseudoscalar and show how relevant the 
Higgs contribution may become. Finally we study the cNMSSM (with universal soft terms at 
the GUT scale), and conclude with a short summary. 

In Fig. 1, we plot the supersymmetric contribution a^ USY as a function of tan /? for various 
slepton (sneutrino/smuon) masses. Here we chose, for simplicity, universal soft slepton masses 
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Figure 1: The SUSY contribution to function of tan/3 for various slepton (sneu- 

trino/smuon) masses. 



Msl at low energy. The gaugino soft terms are assumed to be hierarchical (M3 = 3M2 = 
6M1 = 900 GeV) and fi e ff is chosen such that the lighter chargino x + nas a mass of 175 GeV. 
The experimentally allowed la and 2a regions are indicated as an orange and a yellow band, 
respectively. The full violet, blue or red curves (corresponding to /i e // > 0) and dashed violet, 
blue or red curves (corresponding to fi e ff < 0) include all SUSY contributions. Next to them 
we plot as dot-dashed lines the one- loop chargino/sneutrino contribution separately (corrected 
by large QED logarithms; for sneutrino/smuon masses of 1 TeV they differ hardly from the 
full SUSY contributions). As we mentioned in the previous section, the chargino/sneutrino 
contribution is proportional to the Yukawa coupling of the muon and thus to tan (3. Since this 
contribution obviously dominates the total SUSY contribution (the small difference is mainly 
due to the neutralino diagram), the total SUSY contribution is also roughly proportional to 
tan/3. 

The chargino contribution to carries the same sign as the SUSY parameter /1 (/i e // in 
the NMSSM). Hence, the case fi < (dashed curves in Fig. 1) is disfavoured by the sign of the 
difference of the BNL result and the SM ( 11.91) . When /1 is positive (full curves in Fig. 1), the 
SUSY contribution to is able to account for the 3a deviation. When sleptons are heavy (the 
red curve), however, large values of tan/3 are necessary. On the other hand, for light sleptons 
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(violet curve) the chargino contribution can be large enough even for moderate values of tan /3, 
or even too important when tan /3 is large. 
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Figure 2: Exclusion plot in the tan/3/slepton mass plane for a light chargino of 175 GeV. 

In Fig. 2 we present an exclusion plot in the tan/3/slepton mass plane for a light chargino 
of 175 GeV. In the red and yellow hatched domains the muon magnetic moment differs from 
the BNL result at the 2a and la level, respectively. (In the excluded dark hatched region a 
slepton would be the LSP, or violate constraints from its non-observation at LEP.) Two different 
domains are excluded by a^. 

• In the bottom right-hand corner of Fig. 2, tan/3 is large and sleptons are light. In this 
case, the chargino/sneutrino contribution is strongly enhanced and becomes too large. 

• When sleptons are heavy and tan/3 is low (top and left-hand side of Fig. 2), the chargino 
contribution to is suppressed, and SUSY cannot explain the BNL result. 

In order to see the impact of the chargino mass, we consider the case of a heavier chargino 
of 450 GeV in Fig. 3. Here we chose M3 = 3M2 = 6M\ = 1.5 TeV. One sees that light sleptons 
and large values of tan/3 are required in order to account for the BNL result. 
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Figure 3: Exclusion plot in the tan/3/slepton mass plane for a heavier chargino of 450 GeV. 



These results confirm that the chargino/sneutrino contribution can explain the 3a deviation 
between a^ M and a^ xp in the NMSSM as well as in the MSSM, provided that the supersym- 
metric particles are sufficiently light or tan (3 is large. 

Next we consider the NMSSM with a light pseudoscalar a\ in the spectrum. From the 
formulae (12.211) and (12.241) one finds that its contribution to depends essentially on the 
pseudoscalar mass m ai and its coupling to (down) fermions X d = P' X1 tan (3, where P[ x describes 
the mixing of a\ with the Higgs-doublet sector. Subsequently we chose parameters in the Higgs 
sector such that P' X1 remains approximately constant ~ 0.52; then Xd is entirely determined by 
tan (3. 

In Fig. 4, we plot the contribution to a M which originates from the NMSSM Higgs sector 
only against tan/5 (lower axis) or Xd (upper axis) for m ai = 6.5 GeV; then the contributions 
in section 2.2 are dominated by the light pseudoscalar. We indicate separately the negative 
one-loop contribution (green curve), and the positive two-loop contribution (red curve). For 



m 



6.5 GeV the two- loop contribution dominates, so that the total result (black curve) has 
the same sign as the desired contribution (jl.9p . The contribution behaves as tan 2 /9 (or X%) 
and we find that the Higgs contribution alone can reduce the discrepancy (11.91) below the 2a 
(la) level for X d > 15 (22.5). 

In Fig. 5 we show the contribution to a M from the NMSSM Higgs sector as a function of 
m ai for various values of Xd- As we mentioned in the previous section, the negative one-loop 
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Figure 4: Higgs contribution to function of tan (3 or Xd for m ai = 6.5 GeV. 

contribution dominates below m ai ~ 3 GeV, while the two-loop diagram dominates for larger 
masses. Since both contributions decrease when the pseudoscalar becomes heavy, the total 
Higgs contribution becomes maximal around ~ 6 — 7 GeV. 

Constraints from B physics (essentially from B s — > fi + fi~ and AM S ) depend mainly on the 
loop-induced b — s — a x coupling [30] , and are particularly strong for m ai ~ M Bd s ~ 5 GeV. 
For fixed (and sfermion masses and trilinear couplings at 1 TeV, A = 0.3, fi e ff = 200 GeV, 
and M3 = 3M2 = 6M1 = 1.2 TeV) we varied the remaining parameters tan/5, k and A K in 
the Higgs sector of the NMSSM and obtained regions near m ai ~ 5 GeV, which are always 
excluded - these regions are indicated as dashed parts of the curves in Fig. 5. For larger 
Xd and corresponding larger values of tan the forbidden regions are larger. For X^ = 24, 
the forbidden region extends up to m ai ~ 6.5 GeV. The region m ai > 6.5 GeV, where the 
total Higgs contribution can be relatively large, can be consistent with B physics constraints, 
however. 

The previous results concern certain regions of the relatively large parameter space of the 
general (low energy) NMSSM. Next we consider the impact of a M on the cNMSSM with universal 
scalar and gaugino masses (m and Mi/ 2 ) as well as universal trilinear couplings (A ) at the 
GUT scale. Such a model is motivated by flavour-blind, supergravity-induced SUSY breaking 
scenarios. 

A recent study of the cNMSSM parameter space spanned by A, m , M 1 / 2 and A [39] showed 



12 



5a Hlggs *l0 10 




20 

rn (GeV) 



Figure 5: Higgs contribution to function of the pseudoscalar mass. 

that LEP constraints on Higgs scalars and, notably, a dark matter relic density in agreement 
with WMAP constraints require A < 1, mo <C M\/ 2 and that Aq is determined in terms of M\/2- 
Consequently, the complete Higgs and sparticle spectrum depends essentially only on Mi/ 2 , 
which can a priori vary between 400 GeV and 2-3 TeV (where all other constraints on sparticle 
masses and from B-physics are satisfied; the lower limit on Mi/ 2 originates simultaneously from 
the lower experimental bound on stau masses and LEP constraints on Higgs scalars). 

Using the most recent version of NMSSMTools [35], we have computed the spectrum and 
the supersymmetric contributions to a M in the cNMSSM as a function of M1/2, with the result 
shown in Fig. 6. 

We see that the constraint from a M confines the allowed range of M\/2 to M1/2 < 1 TeV at 
the 2a level, and to 400 GeV < M 1 / 2 < 700 GeV (where the sparticle spectrum is not too 
heavy) at the la level. In fact, the present experimental value could be matched to arbitrarily 
high precision, and a more precise measurement of a M could determine M1/2 completely. 

In the cNMSSM, the leading SUSY contributions to a M originate from the MSSM-like one- 
loop diagrams. Besides the chargino/sneutrino loop, the bino/smuon contribution is also quite 
significant. On the other hand, contributions from the Higgs sector are always negligible in the 
cNMSSM, since the lightest pseudoscalar is not very light. 

To summarize, the deviation ( 11. 9ft of the measured value of a M from the SM can be explained 
in the NMSSM as in the MSSM since - in the absence of a light CP-odd Higgs scalar - the 
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Figure 6: 5a,, function of M 1/2 in the cNMSSM 



corresponding contributions are practically the same. In the presence of a light CP-odd Higgs 
scalar the NMSSM specific contributions to are not negligible in general (depending on Xd 
and m ai , see Figs. 4 and 5) and, for m ai > 3 GeV, can allow for a heavier sparticle spectrum. 
The mSUGRA-like cNMSSM is consistent with constraints from a M , which can help to confine 
the remaining free parameter Mu2- 
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Appendix: Loop functions 

The functions appearing in the one-loop contributions (12.11 12.21) are given in [11]: 

2 

F?(x) = j [l - 6x + 3x 2 + 2x 3 - 6a; 2 In x] 

[I — x) 

F ?(%) = 3 \ q [1 -x 2 + 2xlna:] 
(1 — x) 6 

2 

FF(x) = — [2 + 3a; - Qx 2 + x 3 + Qx In x] 

(1 — x) 4 

F 2 c (x) = - 2(1 ^ )3 [3-4x + x 2 + 21nx], 

The functions appearing in the two-loop contributions (12.131 12.171 2.23, 12.24p are as in [9]: 

, / s z f 1 x(l — x) , x(l — x) , z r „ , / m 

//(,) ^ ~2 / , (1 - g) - z ln cfa = - [2 + In , - /„(«)] 

/ps(^) = z / — r In dx 

J x(1 - x) - z z 



2z 



^T=Iz~ 



Li2 n-I-S-Li 2 fi 1 + ^ 



2z J \ 2z 



/• 1 l-2x(l-a:) 1 xO-x) 

M«) = -2 / — r: r In dx 

J x(l - x) - z z 

= (2z-l)f PS {z)-2z(lnz + 2) 

References 

[1] G . W. Bennett et al [Muon G-2 Collaboration], Phys. Rev. D 73 (2006) 072003 
|arXiv:hep-ex/0602035] . 



[2] F. Jegerlehner, Nucl. Phys. Proc. Suppl. 162 (2006) 22 |arXiv:hep-ph/0608329] . 



[3] K. Hagiwara, A. D. Ma rtin, D. Nomura and T. Teubner, Phys. Lett. B 649 (2007) 173 
|arXiv:hep-ph/0611l"02] . 



[4] M. Davier, Nucl. Phys. Proc. Suppl. 169 (2007) 288 [arXiv:hep-ph/0701163| 



[5] F. Jegerlehner, Acta Phys. Polon. B 38 (2007) 3021 |arXiv:hep-ph70 703125|. 
[6] Z. Zhang, "Muon g-2: a mini review," larXiv:0801.4905l [hep-ph], 

[7] R . Alemany, M. Da vier and A. Hocker, Eur. Phys. J. C 2, 123 (1998) 
|arXiv:hep-ph/9703220] . 



15 



[8] M. Davier, S. Eidelma n, A. Hocker and Z. Zhang, Eur. Phys. J. C 31, 503 (2003) 
|arXiv:hep-ph/0308213] . 



[9] D. Stockinger, J. Phys. G 34 (2007) R45 [arXiv:hep-ph/0 609168|. 

[10] P. Fayet, in Unification of the Fundamental Particles Interactions, eds. S. Ferrara, J. Ellis, 
P. van Nieuwenhuizen, (Plenum, NY, 1980), p. 587; 
J. A. Grifols and A. Mendez, Phys. Rev. D 26 (1982) 1809; 
J. R. Ellis, J. S. Hagelin and D. V. Nanopoulos, Phys. Lett. B 116 (1982) 283; 
R. Barbieri and L. Maiani, Phys. Lett. B 117 (1982) 203; 
D. A. Kosower, L. M. Krauss and N. Sakai, Phys. Lett. B 133 (1983) 305; 
T. C. Yuan, R. Arnowitt, A. H. Chamseddine and P. Nath, Z. Phys. C 26 (1984) 407; 
J. C. Romao, A. Barroso, M. C. Bento and G. C. Branco, Nucl. Phys. B 250 (1985) 295; 
I. Vendramin, Nuovo Cim. A 101 (1989) 731; 

S. A. Abel, W. N. Cottingham and I. B. Whittingham, Phys. Lett. B 259 (1991) 307; 
U. Chattopadhyay and P. Nath, Phys. Rev. D 53 (1996) 1648 |arXiv:hep-ph/9507386] ; 
T. Moroi, Phys. Re v. D 53 (1996) 6565 [Erratum-ibid. D 56 (1997) 4424] 
|arXiv:hep-ph/95 12396] . 

[11] S. P. Martin and J. D. Wells, Phys. Rev. D 64, 035003 (2001) |arX iv:hep-ph/0103067] . 

[12] M. Byrne, C. Kolda and J. E. Lennon, Phys. Rev. D 67 (2003) 075004 
|arXiv:hep-ph/02080"67] . 

[13] S. P. Martin and J. D. Wells, Phys. Rev. D 67 (2003) 015002 [arXiv:hep-ph/0209309| . 

[14] G. Degrassi and G. F. Giudice, Phys. Rev. D 58, 053007 (1998) |[arXiv:hep-ph/9803384] . 

[15] D. Chang, W. F. C hang, C. H. Chou and W. Y. Keung, Phys. Rev. D 63 (2001) 091301 
[arXiv:hep-ph/00092"92] . 

[16] K. M. Cheung, C. H . Chou and O. C. W. Kong, Phys. Rev. D 64, 111301 (2001) 
[arXiv:hep-ph/0103183] . 

[17] C. H. Chen and C. Q. Geng, Phys. Lett. B 511 (2001) 77 |arXiv:hep-ph/0104151| . 

[18] A. Arhrib and S. Baek, Phys. Rev. D 65, 075002 (2002) |arXiv:hep-ph/0104225] . 

[19] S. Heinemeyer, D. S tockinger and G. Weiglein, Nucl. Phys. B 690, 62 (2004) 
|arXiv:hep-ph/03 12264] . 

[20] S. Heinemeyer, D. St ockinger and G. Weiglein, Nucl. Phys. B 699, 103 (2004) 
|arXiv:hep-ph/0405255] . 

[21] T. F. Feng, X. Q. Li, L Lin, J. Maalampi and H. S. Song, Phys. Rev. D 73 (2006) 116001 
|arXiv:hep-ph/0604l"7l] . 



16 



[22] T. F. Feng, L. Sun and X. Y. Yang, "Electroweak and supersymmetric two-loop corrections 
to lepton anomalous magnetic and electric dipole moments," larXiv:0805.1122l [hep-ph]. 

[23] J. F. Gunion, D. H ooper and B. McElrath, Phys. Rev. D 73, 015011 (2006) 
|arXiv:hep-ph/0509024] . 

[24] V. Barger, C. Kao, P. Langacker and H. S. Lee, Phys. Lett. B 614 (2005) 67 
|arXiv:hep-ph/0412136] . 

[25] H.P Miles, M. Srednicki and D. Wyler, Phys. Lett. B 120 (1983) 346; 
J.M. Frere, D.R. Jones and S. Raby, Nucl. Phys. B 222 (1983) 11; 
J.P Derendinger and C.A. Savoy, Nucl. Phys. B 237 (1984) 307; 

J.R. Ellis, J.F. Gunion, H.E. Haber, L. Roszkowski and F. Zwirner, Phys. Rev. D 39 (1989) 
844; 

M. Drees, Int. J. Mod. Phys. A 4 (1989) 3635. 

[26] B. A. Dobrescu, G. L andsberg and K. T. Matchev, Phys. Rev. D 63 (2001) 075003 
|arXiv:hep-ph/0005308| ; 

B. A. Dobrescu and K. T. Matchev, JHEP 0009 (2000) 031 [arXiv: hep-ph /00081 92| . 



[27] U. Ellwanger, J. F. Gunion, C. Hugonie and S. Moretti, |arXiv: hep-ph/0305 1 09] (in 
"Physics interplay of the LHC a nd the ILC", G. Weigl ein et al. [LHC/LC Study Group], 
Phys. Rept. 426 (2006) 47); |arXiv:hep-ph/0401228| (in "The Higgs working group: 
Summary report 2003", K. A. Assamagan et al. [Higgs Working Group Collaboration], 
arXiv:hep-ph/0406152|) ; 



R. Dermisek and J. F. Gunion, P hys. Rev. Lett. 95 (2005 ) 041801 |arXiv:hep-ph/0502105| , 
Phys. Rev. D 73 (2006 ) 111701 [arXiv:hep-ph/05 10322] , Phys. Rev. D 75 (2007) 075019 
|arXiv:hep-ph/0611142| ; 

U. Ellwanger, J. F . Gunion and C. Hugonie, JHEP 0507 (2005) 041 
|arXiv:hep-ph/0503203] ; 

S. Chang, P. J. Fox and N. Weiner, JHEP 0608 (2006) 068 [arXiv:hep-ph/0511250] ; 

P. W. Graham, A. Pierce and J. G. Wacker, "Four taus at the Tevatron," 



arXiv:hep-ph/0605162 



S. Moretti, S. Munir and P. Poulose, Phys. Lett. B 644 (2007) 241 |arXiv:hep-ph/0608233| ; 
S. Chang, P. J. Fo x and N. Weiner, Phys. Rev. Lett. 98 (2007) 111802 
|arXiv:hep-ph/0608310] ; 

T. Stelzer, S. Wiesen feldt and S. Willenbrock, Phys. Rev. D 75 (2007) 077701 
|arXiv:hep-ph/0611242] ; 



U. Aglietti et al, "Tevatron-for-LHC report: Higgs," |arXiv:hep-ph/0612172[ 

E. Fullana and M. A. Sanchis-Lozano, Phys. Lett. B 653 (2007) 67 |arXiv:hep-ph/0702190| ; 
K. Cheung, J. Song and Q. S. Yan, Phys. Rev. Lett. 99 (2007) 031801 
|arXiv:hep-ph/0703149] ; 

M. A. Sanchis-Lozano, "A light non-standard Higgs boson: to be or not to be at a (Super) 
B factory?," larXiv:0709.3647l [hep-ph]; 

M. Carena, T. Han, G. Y. Huang and C. E. M. Wagner, JHEP 0804 (2008) 092 



17 



|arXiv:0712.2466l [hep-ph]]; 

J.R. Forshaw, J.F. Gunion, L. Hodgkinson, A. Papaefstathiou and A.D. Pilkington, 
larXiv:0712.35T0l [hep-ph]; 

Z. Heng, R. J. Oakes, W. Wang, Z. Xiong and J. M. Yang, "B meson Dileptonic Decays 

in NMSSM with a Light CP-odd Higgs Boson," larXiv: 0801 .11691 [hep-ph]; 

A. Djouadi et al, "Benchmark scenarios for the NMSSM". larXiv: 0801 .43211 [hep-ph]; 

X. G. He, J. Tandean and G. Valencia, "Rare Decays with a Light CP-Odd Higgs Boson 

in the NMSSM," larXiv:0803.4330l [hep-ph]. 

[28] M. Krawczyk, Acta Phys. Polon. B 33, 2621 (2002) | arXiv:hep-ph/0208076| . 

[29] S. Schael et al. [ALEPH, DELPHI, L3 and OPAL Collaborations], Eur. Phys. J. C 47 
(2006) 547 

[30] F. Domingo and U. Ellwanger, JHEP 0712 (2007) 090 [arXiv:0710.3"7T4l [hep-ph]]. 

[31] S. Laporta and E. Remiddi, Phys. Lett. B 379, 283 (1996) |arXiv:hep-ph/9602417 | , Phys. 
Lett. B 301, 440 (1993); 

T. Aoyama, M. Hayakawa, T. Kinoshita and M. Nio, Phys. Rev. Lett. 99, 110406 (2007) 
|arXiv:0706.3l96l [hep-ph] ] ; 

T. Kinoshita and M. Nio, Phys. Rev. D 70, 113001 (2004) |arXiv:hep-ph/0402206| . 

[32] T. Kinoshita and M. Nio, Phys. Rev. D 73, 053007 (2006) |arXiv:hep-ph/0512330| . 

[33] J. Bijnens and J. Prades, Mod. Phys. Lett. A 22 (2007) 767 |arXiv:hep-ph/0702170| . 

[34] A. Czarnecki, W. J. Marcia no and A. Vainshtein, P hys. Rev. D 67 (2003) 073006 [Erratum- 
ibid. D 73 (2006) 119901] |arXiv:hep-ph/0212229 . 



U. Ellwanger, J. F . Gunion and C. Hugonie, JHEP 0502 (2005) 066 
|arXiv:hep-ph/0406215] ; 

U. Ellwanger and C. Hugonie, Comput. Phys. Commun. 177 (2007) 399, 
|http://www.th.u-psud.fr/NMHDECAY/nmssmtools.html 



[36] P. Skands et a/., JHEP 0407 (2004) 036 |arXiv:hep-ph/0311123 



[37] B. Allanach et al, "SUSY Les Houches Accord 2." larXiv:0801.0045l [hep-ph]. 

[38] J. P. Leveille, Nucl. Phys. B 137 (1978) 63. 

A. Djouadi, U. Ellwanger and A. M. Teixeira, "The constrained next-to-minimal super- 
symmetric standard model," larXiv :0803.0253 [hep-ph]. 



18 



